STABILITY OF VORTEX SHEETS AND ENTROPY WAVES 
IN STEADY SUPERSONIC EULER FLOWS OVER LIPSCHITZ WALLS 
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Abstract. We establish the weU-posedness of compressible vortex sheets and entropy waves in two- 
dimensional steady supersonic Euler flows over Lipschitz walls under a BV boundary perturbation. 
In particular, when the total variation of the incoming flow perturbation around the background 
strong vortex sheet /entropy wave is small, we prove that the two-dimensional steady supersonic 
Euler flows containing a strong vortex sheet/entropy wave past a Lipschitz wall are L^-stable. 
Both the Lipschitz wall (whose boundary slope function has small total variation) and incoming 
flow perturb the background strong vortex sheet /entropy wave. The weak waves are reflected after 
nonlinear waves interact with the strong vortex sheet/entropy wave and the wall boundary. Using 
the wave-front tracking method, the existence of solutions in BV over Lipschitz walls is first shown, 
when the incoming flow perturbation of the background strong vortex sheet/entropy wave has small 
total variation. Then we establish the -contraction of the solutions with respect to the incoming 
flows. To achieve this, a Lyapunov functional, equivalent to the L'^-distance between two solutions 
containing strong vortex sheets/entropy waves, is carefully constructed to include the nonlinear 
waves generated both by the wall boundary and from the incoming flow. This functional is then 
shown to decrease in the flow direction, leading to the L^-stability, as well as the uniqueness, of the 
solutions. Furthermore, the uniqueness of solutions extends to a larger class of viscosity solutions. 



1. Introduction 

We study the well-posedness of two-dimensional steady supersonic Euler flows past a curved Lips- 
chitz wall containing strong vortex sheets/entropy waves in the L^-norm. The inviscid compressible 
flows are governed by the two-dimensional steady Euler system: 

'ipu)x + {pv)y = 0, 
{pu^ +p).j, + {puv)y = {), 
{puv)x + {pv"^ + p)y = 0, 



t(pu(ii; + 2)),. + (H^ + 2)), = 0, 

with {u,v), p, p, and E representing the fluid velocity, scalar pressure, density, and total energy, 
respectively. Furthermore, the total energy E is explicitly given by 

E=^{u^+v^)+e{p,p), 

where the internal energy e can be written as a function of {p,p) defined through the thermody- 
namical relations. The temperature T and entropy S are the other two thermodynamic variables. 
In the case of an ideal gas, the pressure p and internal energy e can be expressed as 

p = RpT, e = CyT (1.2) 

with the adiabatic index 7 given by 

7 = 1 + — >1. (1.3) 
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In particular, in terms of the density p and entropy S", we have 

p = p{p^ S) = Kp^e^/'% e = -^p^-^e^l"'' = (1.4) 

7 — 1 7 ~ 1 

The constants R,c^, and k in the above relations are all greater than zero. 

When the entropy S = constant, the flow is called isentropic. In this case, the pressure p can be 
written as a function of the density p, p = p{p), and the flow is governed by the isentropic Euler 
equations: 

{pu)x + {pv)y = 0, 
< {pu"^ +p)r, + {puv)y = 0, (1.5) 

^ {puv):r + (pw^ + p)y = 0. 

Then, by scaling, the pressure-density relation is 

p{p) = —■ (1-6) 

7 

The adiabatic exponent 7 > 1 corresponds to the isentropic polytropic gas. The limiting case 7 = 1 
corresponds to the isothermal flow. Define 

c= \jpp{p, S) 

as the sonic speed. For polytropic gases, the sonic speed is c = yyp/p- 

The flow type is classified by the Mach number M = v^"^+"^ , When M > 1, system (jl.ip or 
()1.5p governs a supersonic flow (i.e., + v'^ > c^), which has all real eigenvalues and is hyperbolic. 
For M < 1, system (jl.ip or (jl.Sp governs a subsonic flow (i.e., u'^ + v'^ < c^), which has complex 
eigenvalues and is elliptic- hyperbolic mixed and composite. When M = 1, the flow is called sonic. 

We are interested in whether compressible vortex sheets/entropy waves in supersonic flow over the 
Lipschitz wall are always stable under the BV perturbation of the incoming flow. Multidimensional 
steady supersonic Euler flows are important in many physical applications (cf. Courant-Friedrichs 
|11|). In particular, when the upstream flow is a uniform steady flow above the plane wall in 
x < all the time, the flow downstream above a Lipschitz wall in x > is governed by a steady 
Euler flow after a sufficiently long time. Moreover, compressible vortex sheets and entropy waves 
occur ubiquitously in nature and are fundamental waves. Furthermore, since steady Euler flows 
are asymptotic states and may be global attractors of the corresponding unsteady Euler flows, it is 
important to establish the existence of steady Euler flows and understand their qualitative behavior 
to shed light on the long-time asymptotic behavior of the unsteady compressible Euler flows, one 
of the most fundamental problems in mathematical fluid dynamics which is still wide open. 

We observe that the stability of contact discontinuities for the Cauchy problem for strictly hyper- 
bolic systems in one space dimension under a BV perturbation has been studied by Sable- Tougeron 
|24j and Corli-Sable- Tougeron [12]. In particular, the reflection coefficients, such as Kn here, are 
required to be less than one, which is the stability condition for the mixed problem in the strip 
{{t,x) : t > 0, — 1 < X < 1} in the earlier works; see, e.g.. Sable- Tougeron [2l]. Working with the 
non-isentropic Euler system (jl.ip and a uniform upstream flow, Chen-Zhang-Zhu |10] first proved 
the global existence in BV of supersonic Euler flows containing a strong vortex sheet /entropy wave 
under the BV perturbation of the Lipschitz wall by using the Glimm scheme. The essential dif- 
ference between system (jl.ip as analyzed in [10] (and in Sections 2-7 here) and strictly hyperbolic 
systems as considered in |12j is that two of the four characteristic eigenvalues coincide and have 
two linearly independent eigenvectors which determine precisely the compressible vortex sheets and 
entropy waves so that two independent parameters are required to describe them, respectively. 

In this paper, for completeness, we first show, via the wave-front tracking method, the existence 
of solutions to the problem when a small BV perturbation is added to the uniform incoming 
flow. Then the L^-stability of entropy solutions containing strong vortex sheets/entropy waves is 
established. As corollaries of these results, the estimates on the uniformly Lipschitz semigroup S 
of entropy solutions generated by the wave-front tracking approximations are obtained, and the 
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uniqueness of weak solutions containing strong vortex sheets/entropy waves is established in a 
larger set of solutions, namely the class of viscosity solutions. 

In the following, we focus mainly on the problem in the region over the Lipschitz wall for the 
supersonic Euler flows U (u, v,p, p) governed by system (jl.ip . given that the corresponding problem 
for the isentropic system (jl.Sp is simpler to analyze. The subsequent figure provides the schematic 
diagram for the problem we study: 



Vortex sheet/Entropy wave 




y=g(x) 



Figure 1.1. Stability of the compressible vortex sheet /entropy wave in supersonic flow 

The boundary and initial data in the problem are as follows: 

(i) There is a Lipschitz function g £ Lip(]R+;]R) such that 

£,(0) = g'{0+) = 0, lim arctan(c/'(a;+)) =0, g' £ BV{R+;M) 

and 

TV{g'{-)) < e for some constant e > 0. 
Denote := {{x^y) : y > g{x),x > 0}, F := {{x,y) : y = g{x),x > 0}, and n(x±) = 
( 9 ^]^g outer normal vectors to F at the respective points xit (cf. Fig. 11.11) . 

V(9'(^±))2 + l _ ^ f f K & V 

(ii) The incoming flow U = U{y) := Uq + Uq aA x = {) is composed of two parts: 

(a) The upstream flow Uq consists of one straight vortex sheet /entropy wave y = y^ > 



and two constant vectors U, 



satisfying 







U- , when < y < y^, and Uq 



[/+, wheny > y,^ > 0, 



0, 



u± > c± > 0, 



where c± = \rw±/p± is the sonic speed of state U±. 
(b) The perturbation Uq = iuo,vo,Po, po)iy) G n By(M;M^) at x = so that 
TV([/o) < 1. 

Then we consider the following initial-boundary value problem for system (jl.ip : 

Cauchy Condition : U\^=o = U{y) = + Uq\ (1.7) 

Boundary Condition : (n, • n = on F. (1-8) 

Definition 1.1 [Admissible entropy solutions). A BV function U = U{x,y) is said to be an 
entropy solution of the initial-boundary value problem (jl.ip and ()1.7p - ()1.8p if and only if the 
following conditions hold: 

(i) [/ is a weak solution of (2.1) and satisfies 

U\x=o = U{y) and {u,v) ■ n|y=g(^) = in the trace sense; 

(ii) U satisfies the steady entropy Clausius inequality: 

{puS)x + {pvS)y > (1.9) 
in the distributional sense in il. including the Lipschitz wall boundary. 
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One of the essential developments within this paper is to develop suitable methods to deal with 
the challenges caused by the nonstrictly hyperbolicity of the system and the Lipschitz wall boundary, 
in comparison with the previous progress with the strictly hyperbolic systems of conservation laws, 
particularly to the analysis of the Cauchy problem. For supersonic Euler flow with a strong shock- 
front emanating from the wedge vertex, Chen-Li worked out the issue for a Lipschitz wedge 
boundary. We now discuss some main differences in our work here from the Cauchy problem and 
the resulting key difficulties. We remark that, in the case of the Cauchy problem concerning only 
weak waves, the decrease of the Lyapunov functional and the L^-stability of the solutions were 
obtained through the cancellation of distances on both sides of waves. In the presence of a strong 
shock, for the L^-stability of solutions of the Cauchy problem for strictly hyperbolic systems of 
conservation laws, the Lyapunov functional was found to decrease by employing the strength of the 
strong shock to control the strengths of weak waves of the other families (e.g., see Lewicka-Trivisa 
|21j). In contrast with our Lipschitz wall problem, which is a problem of initial-boundary value 
type, there is no such cancellation by the boundary as only one-side is possible near it. Furthermore, 
no strong vortex sheets/entropy waves (characteristic discontinuities) nor strong shocks are present 
to handle the strength of the weak waves of the other families, and the terms in the estimates 
for the first and fourth family carry different signs. As such, it is difficult to say whether the 
functional can be made to decrease for our case of strong vortex sheets and entropy waves with 
multiplicity of eigenvalues. One of the key steps to resolve this is to use the physical feature of 
the boundary condition that the flow of two solutions near the boundary must run in parallel (also 
see [S]). This observation helps us to obtain additional quantitative relations near the boundary. 
Then, applying suitable weights and adjustments in the coefficients of the Lyapunov functional and 
using the cancellation between the different families, the functional is found to decrease in the flow 
direction. 

The rest of the paper is organized as follows. In Section 2, we recall some fundamental properties 
of the two-dimensional steady Euler system (jl.ip and discuss related nonlinear waves and wave 
interaction estimates. In Section 3, the wave- front tracking algorithm is discussed, working in the 
presence of strong vortex sheets/entropy waves, the suitable interaction potential Q is constructed, 
including the effect of the Lipschitz wall, and the existence of entropy solutions in BV is established 
for the initial-boundary value problem. In Section 4, we construct the Lyapunov functional # 
(equivalent to the L^-distance between two entropy solutions U and V) to include the nonlinear 
waves produced by the wall boundary vertices. Then, in Section 5, the monotone decrease of 
the functional ^ is established in the flow direction, leading to the L^~stability of the solutions 
containing strong vortex sheets/entropy waves. Using the estimates established in Sections 3-5, 
in Section 6, we obtain the existence of a Lipschitz semigroup of solutions generated by a wave- 
front tracking approximation, as well as some estimates on the uniformly Lipschitz semigroup S 
produced by the limit of wave- front tracking approximations. Moreover, the uniqueness of solutions 
with strong vortex sheets/entropy waves is obtained in the larger class of viscosity solutions. 

2. Adiabatic Euler equations: Nonlinear waves and wave interactions 

In this section, we first present some basic properties of the steady Euler system (jl.ip . Then 
related nonlinear waves and interaction estimates are discussed, which will be employed in the later 
sections. 

Consider the following vector functions of the solution U: 

W{U) = {pu,pu^ +p,puv,pu{h H ))'^, 

H{U) = {pv, puv, pv^ + p, pv{h H 

where h = ^-^Z^^p ■ Then the steady Euler equations in (jl.ip can be expressed in the following 
conservative form: 

W{U), + H{U)y = 0, U = {u,v,p,p)'^ (2.1) 



Li-STABILITY OF VORTEX SHEETS AND ENTROPY WAVES 5 

When U{x,y) is a smooth solution, system (j2.ip is equivalent to 

VuW{U)U, + VuH{U)Uy = 0. (2.2) 
Then the roots of the fourth degree polynomial 

det{XVuW{U) - VuH{U)), (2.3) 
are the eigenvalues of ()2.ip : that is, the solutions of the equation 

{v - Xuf{{v - Xuf - c^{l + A^)) = 0, (2.4) 

where c = is the sonic speed. For supersonic flows (i.e. + v"^ > c^), system ()2.ip is 

hyperbolic. Specifically, when u > c, system (|2.ip has four real eigenvalues in the x-direction: 

UV + {—l)'^cVu'^ + 1)2 — 

Xd = 9 9 , "=1,4; 

Xk = -, k = 2,3, (2.5) 
u 

with the four corresponding linearly independent eigenvectors given by 

rd = Kd{-Xd,l,p{XdU-v), 2 ) , a = 1,4, 

r2 = (n,^;,0,0)^, = (0, 0, 0, p)^, (2.6) 

where Kd the re-normalization factors such that r^^- VA^ = 1, given that the dth-characteristic fields, 
d = 1,4, are genuinely nonlinear. The second and third linearly degenerate characteristic fields 
satisfy ■ VA^ = 0, k = 2,3, which correspond to vortex sheets and entropy waves, respectively. 

The wave curves in the phase space are now described. The Rankine-Hugoniot jump conditions 
for (QA^ are 

a[Wiu)] = [Hiu)], (2.7) 

and the discontinuity propagates with the speed a. 

There are two different waves associated with the fields Afc = fc = 2, 3, with the corresponding 
linearly independent right eigenvectors r^jk = 2,3, in (|2.6p : 

Vortex sheets: 



Entropy waves: 



C2{Uo): a = - = ^, p = po, S = Sq, + ^ 4 + vl (2.8) 

U Uq 

CsiUo) : a = - = —, p = po, {u,v) = {uo,vo), 5 / Sq. (2.9) 

U Uq 

Albeit the two contact discontinuities, the vortex sheet and the entropy wave, above match as a 
single discontinuity in the physical (x, y)-plane, two independent parameters are needed to describe 
them in the phase space U = {u,v,p,p) since there are two linearly independent eigenvectors 
corresponding to the repeated eigenvalues A2 = A3 = ^ of the linearly degenerate characteristics 
fields. 

The nonlinear waves associated with Xd,d = 1,4, are shock waves and rarefaction waves. The 
shock waves have their speeds of propagation given by 



UqVq + {-lYco^Jul +vl- Cn 



where Cn = ^^-7- and 60 = — "V^TT- Substituting ad, d = l,A, into (j2.7p . the d-Hugoniot curve 



73. 

a — ad-— 5 lo , a — i,4 

— Cq 

_ £0 

Sd{Uo) through the state Uq is 

Sd{Uo) : [p] = 7^[/o], [u] = -ad[v], poiadUo - Vo)[v] = [p], d = 1,4. 

Written as ^^"(C/o), d = 1,4, the half curves of Sd{Uo) for p > po in the phase space are said to 
be the shock curves on which any state forms a shock with the below state Uq in the {x, y)~plane 
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respecting the entropy condition (jl.9p . Furthermore, for each d = 1 or d = A, the curves S'^{Uo) 
and R'^{Uq) at the state Uq have the same curvature. 

If [/ is a piecewise smooth solution (see also [TO]), then any of the following conditions below is 
equivalent to the entropy inequality ()1.9p in Definition 1.1 for a shock wave: 

(i) The physical entropy condition: The density increases across the shock in the flow direction, 

Pback < Pfront- (2-10) 

(ii) The Lax entropy condition: On the dth-shock, the shock speed ad satisfies 

A(i(back) <ad< Ad(front) for d = 1,4, (2.11) 
CTi < A2,3(back), A2,3(front) < 0-4. (2.12) 

The rarefaction wave curves R^{Uq) through the state Uq in the state space are given by 

Rd '■ dp = c^dp, du = —Xddv, p{XdU — v)dv = dp for p < po) d = l,4. (2.13) 

We next discuss several essential properties of the nonlinear waves and related wave interaction 
estimates in Lemmas 2.1-2.7 below. These facts will be used in the later sections. We also refer 
the reader to Chen-Zhang-Zhu [TO| for further details. 

2.1. Riemann Problems and Riemann Solutions 

We focus on the related Riemann problems and their solutions in this section, which serve as the 
building blocks for the front tracking algorithm for the initial-boundary value problem (|2.1|) and 

dLZD-dll]). 

Riemann problem of lateral-type. We note that the straight-sided wall problem is the case when 
problem ()2.ip and (|1.7p ~ (|1.8p is considered with the boundary g = 0. It can be seen that, when the 
angle between the straight-sided wall and the flow direction of the incoming flow is zero, problem 
(|2.ip and ()1.7P " ()1.8p has an entropy solution made up of two constants states C/_ = (ii_, 0,p_, /9_) 
and U+ = (ti+, 0,p+, p+), satisfying u± > c± > in the subdomains il+ and of O separated 
by a straight vortex sheet /entropy wave. These are precisely the states U- and U+ below and 
above the large vortex sheet /entropy wave. The principal aim of this paper is to establish the 
L^-well-posedness for problem (|2.ip and ()1.7p - ()1.8p for the solutions near the background solution 
containing a strong vortex sheet /entropy wave {f/-, with g = 0. 

It has been observed in [TT] that, if the angle between the flow direction of the front state and 
the wall at a boundary vertex is smaller than vr and larger than the extreme angle determined by 
the incoming flow state and 7 > 1, then a unique 4-shock is generated, separating the front-state 
from the supersonic back-state. If the angle between the flow direction of the front-state and the 
wall at a boundary vertex is larger than vr and less than the extreme angle, then a 4-rarefaction 
wave is produced, emanating from the vertex. These waves are easily seen through the shock polar 
analysis (cf. [10\ lllj). This signifies that, when the angle between the flow direction of the front- 
state and the wall at a boundary vertex is close to vr, the lateral Riemann problem can be uniquely 
solved. For further details, see Lemma 2.3 and [10] , For an indepth discussion, we also refer to 
Courant-Friedrichs 1111. 



Riemann problem involving only weak waves. Consider the subsequent initial value problem with 
piecewise constant initial data: 



(2.14) 



with the constant states Ua and Uf, denoting the above state and below state with respect to the 
line y = yo, respectively. Then there is e > so that, for any states Ub, Ua in the neighborhood 
Oe{U+) of [/+, or Ub, Ua in the neighborhood 0^(11^) of ?7_, the initial value problem ()2.14p has a 
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unique admissible solution consisting of four waves, consisting of shocks, rarefaction waves, vortex 
sheets and/or entropy waves. 

Riemann problem involving the strong vortex sheets/entropy waves. From now on, the notation 
{Ub,Ua} = (qi, 04) will be used to write Ua = $ (0:4, 03, 02, ai; J/?,) as the solution of the 

Riemann problem, where $ G C^, and aj is the strength of the j-wave. For any wave with 
Ub G Os{U-) and Ua G 0^{U^), we also use {Uh,Ua} = (0, a"2, era, 0) to denote the strong vortex 
sheet/entropy wave that connects Ub and Ua with strength ((T2,(T3). That is, 

Um. = ^{(72] Ub) ■■= {ube"'^,Vbe"'^,Pb,Pb) , Ua = ^{(73; Um) ■= iUm,Vm,Pm,Pme''^) ■ 

Particularly, we observe that 

U+ = iu+,0,p+,p+) = (tx_e'^2",0,p_,p_e^^°). 

We write G {as, 0-2; Ub) = $3(0-3; ^'2(o"2; Ub)) for any Ub G Oe{U-). Then we have 
Lemma 2.1. The vector function G(as,a2;Ub) satisfies 

Ga, i'J3,<T2;Ub) = (nbe"^7;be"^0,0) , G^, (a3,a2;f/fe) = (0, 0, 0, p^e'^^^), (2.15) 

and 

Vc/G(a3,a2;[/fe) =diag(e^^e'^^l,e^^). (2.16) 

Furthermore, for the plane vortex sheet and entropy wave with the lower state U- = (u_, 0,p_, p_), 
upper state f/+ = (n+, 0,^+, /3+), and strength ((T2o,<73o), 

det{u{U+),G^,{a3o,a2o;U.),G^,{a3o,a2o;U.),VuG{a3o,a2o;U.)-ri{U^)) > 0. (2.17) 

These can be easily obtained from direct calculations and are thus omitted. 

The properties in (|2.15|) - (|2.17|) above play a fundamental role in achieving the necessary estimates 
on the strengths of reflected weak waves in the interaction between the strong vortex sheet / entropy 
wave and weak waves (see the proofs for Lemmas 2.4-2.7). 

2.2. Wave Interactions and Reflection Estimates. In the following, several essential estimates are 
provided on wave interactions and reflections. For their proofs and all the related details, we refer 
to [lU]. 

Weak wave interactions estimates. For the weak wave interaction away from both the strong 
vortex sheet /entropy wave and the wall boundary in the regions or we have the following 
estimate: 

Lemma 2.2. Assume that 

Ub,Um.,UaeOe{U+), Or Ub, U^, Ua € Oe{U^) , 

are three states with {Ub,Um} = (01,02,03,04), {Um.,Ua} = (/?!, /32, /^s, /34)- Then {Ub,Ua} = 
(71,72,73,74) with 

7i = ai + A + 0(l)A(a,/3), (2.18) 
where A(a,/3) = (|a4| + |a3| + |a2|)|/3i| + |a4|(|/32| + |^3|) + I]j=i,4 Aj(a, /?) and 

A>,/)) = |"'„ , "^^"•^'^"^ (Z19) 

l|aj||/3j|, otherwise. 

Estimates on the boundary perturbation of weak waves and the reflection of weak waves on the 

boundary. We write {Gi{ai,bi)}'^Q for the points {(«/, &/)}/^o ™ y)^plane with < ai < a/+i. 
Define 

= arctan , 91 = 61,1+1-61.1,1, 0-1,0 = 0, 

Qi+i = {{x, y) :x e [a/, a^+i] ,y > bi + {x- ai)tan{6 1,1+1)}, (2.20) 
r,+i = {{x,y) : x G {ai,ai+i),y = bi + {x - a/)tan(6l;,i+i)}. 
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and the outer normal vector to F/: 

{k+i - bi,ai - ai+i) I ■ la a \ la a \\ 01 \ 

V(o«+i - o«) + W+i - aiY 
With the constant state C/, consider the following initial-boundary value problem: 

'(2.1) in^^i+i, 

U\x=a, = U, (2.22) 
^{u,v) ■ n;+i =0 on r^+i. 

Lemma 2.3. Suppose {Um,Ua} = (/3i, /32, /^a, 0) and {Ui,Um} = (0,0,0,04) with 

{ui,vi) • n; = 0. 

Then there exists a unique solution f/^+i of problem (|2.22|) such that {Ui^i,Ua} = (0,0,0,^4) and 
Ui+i ■ (n;+i,0, 0) = 0. Moreover, 

(^4 = 04 + KbiPi + Kb2h + Kb3l33 + Kboei , (2.23) 

where Ki,i, Kf,2, -f^bs, and K^q are C'^ -functions of Ps, /32, 04, Oi+i, and Ua satisfying 

-K'bl|{6»j=a4=/3i=/32=/33=0,(7a=C/-} = 1) ^M|{6l;=a4=/3i=/32=/33=0,C/a=C/-} =0' ^ = 2, 3, (2.24) 

and Kbo is bounded. In particular, K^o < at the origin. 

This lemma has two purposes. The first is to estimate the weak waves generated by the vertices on 
the Lipschitz wall boundary. This boundedness will be used to control the boundary perturbation; 
see ()3.2p in the construction of the wave interaction potential Q{x). The second is to estimate 
the strength of the reflected wave ^4 with respect to the incident wave oi. Property ()2.24p of the 
coefficients will play an important role to control the reflected waves. 

Estimates on the interaction between the strong vortex sheet/entropy wave and weak waves from 
below. Estimate ()2.25p below plays a key role in ensuring the L^-stability of entropy solutions, 
especially for the existence of the constants w\ and ^4 in Lemma 5.1 (see below). This estimate 
also ensures the existence of K* E {Ku, 1) in the construction of the wave interaction potential 
Q{x) in (1321). 

Lemma 2.4. Let Ub, Um G Oe{U^) and Ua G Oe{U+) with 

{Ub.Um} = (0,02,03,04), {Um.Ua} = (/?! , Cr2 , ds, 0) . 

Then there exists a unique (5i, cjg, (73, ^4) such that the Riemann problem ()2.14p admits an admissible 
solution that consists of a weak 1—wave of strength 5i, a strong vortex sheet of strength a2, a strong 
entropy wave of strength CJ3, and a weak A— wave of strength 5/^: 

{Ub,Ua} = {5i,a'2,a'^,5i). 



Moreover, 



5i = l3i+ Kiiai + 0(1)A', 5^ = K1404 + 0(1)A', 
4 = 0-2 + 02 + i^i2a4 + 0(1)A', 4 = 0-3 + 03 + ^1304 + 0(1)A', 

A4(f/+)e2'^2o+<T3o - A4(f/_) 



I ^11 1 



{a4=Q!3 ="2=0, (T2='^20iO-3 =0-30} 



A4(C/+)e2'^2o+<^3o + A4(f/- 



< 1, (2.25) 



where A' = |/3i|(|o2| + |a3|), and Ylij=2 l-^iil i^ bounded. 



Lemma 2.5. The coefficient \K\a\( „ / / \ in the strenqth 5a of a weak J-wave 

■/•J I J^'*l|a4=a3=a2=0,tT2=cr20, 0-3=0-30 J- a 4 J t 

in Lemma 2.4 remains bounded away from zero. 

Proof. By Lemma 2.4, we can find a unique solution {6i,a2,cr'^, 64) as a C^-function of 02, 03, 04, /3i, 
(72, ^3, and Ub to 

$4((^4; G{a'^, 4; $i(5; Ub))) = G(ct2, ^3; $i(/3i, $(04, 03, 02, 0; Ub))). (2.26) 
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That is, 

= ^i(«2,a3,a4,/3i,o-2,o-3), i = 2,3; (5j = (5j (a2, as, "4, /?!, 0-2, 0-3), j = l,4, 

where we have omitted Ub for simphcity. 
From |10] . we know that 

1 

Kij = j dajj{a2,a3,9a4^,(3i,a2,cr3)de, j = l,4. 


Differentiate ()2.26p with respect to 04, and let /3i = = 03 = 02 = 0, (T2 = o"2o, and CJ3 = (T30. 
We obtain 

Vc/G(cr30,(T2o;f/-) • r4(C/„) = ^4 r4 + 0-3^^3 ((730, c72o; t/-) 

+ 9a4<T2 G<x2 ((730, 0-20; J/-) + 9^4^! Vc/G(o-3o, o-2o; U-) ■ ri(C/_). 

By Lemma 2.1, we have 

_ det ( V[/ G(cr3o , 0-20 ; C/- ) • r4 ( [/_ ) , G^^ {(T30 , 0-20 ; C^- ), G<^2 (^30 , 0-20 ; f^- ), V^/ G (0-30 , (T20 ; f/- )• ri ( C/_ ) ) 



det(r4 ([/+), G<^3 (0-30, 0-20; f/-), (^^30, 0-20; U^), V[/G(cr3o, CJ20; U-) ■ ri ([/„)) 

Kl(f/_)K4(t/-)/0-'u2 e2<^20+'T30(A4([/_) _ Ai(C/_)) 



Ki(C/_)K4(C/+)p2 ^2 e'^2o+<x3o(A4([/+)e2'^2o+<T3o _ Ai(C/_)) 
2K4(C^-)e'"20A4([/_) 



K4([/+)(A4(C/+)e2'^2o+^3o + \^{U.)) 
This completes the proof. 



> 0. 



Estimates on the interaction between the strong vortex sheet/entropy wave and weak waves from 
above. We have 

Lemma 2.6. Let Ub G Os{U^) and Um, Ua G Os{U+) with 

{Ub,Un,} = (0,^2,^73,04), {Un„Ua} = {^1,^2,^3,0). 

Then there exists a unique 0"2, 1T3, 64) such that the Riemann problem ()2.14p admits an admissible 
solution that consists of a weak 1 — wave of strength 5i, a strong vortex sheet of strength 02, a strong 
entropy wave of strength 03 , and a weak 4— wave of strength 84 : 

{Ub, Ua} = 0-2, 0-3,54). 

Moreover, 

61 = K2i(3i + 0(1) a", 4 = ^2 + /32 + K22P1 + 0(1) a", 

(73 = <T3 + f33 + K23P1 + 0(1)A", 54 = 04 + i^24A + 0(1) A", 

where Yl^j=i l-^2j| is bounded and A" = |a4|(|/32| + l/^sD- 



The constant K21 here is used in the definition of weighted strength of weak waves in ()3.ip . 
Lemma 2.7. The coefficient \K2\\ia a a n 1 / \ in the strenqth 6-1 of a weak 1-wave in 

■IJ I ^J^I|/i3=/i2=/5l=U, 0-2=0-20, cr3=o-3o| ^ ^ ■> 

Lemma 2.6 remains bounded away from zero, while the reflection coefficient |/^24| < 1- 

Proof. By Lemma 2.6, we can find a unique solution (5i, (T2, (T3, J4) as a G2-function of 02,03,04, 
<72, 0-3, and Ub to 

$4(<54; G(a^, 4; $i(<5; Ub))) = m /Ss, /92, A; $4(04; 0(^3, ^2; Ub))). (2.27) 

That is, 

(^i = cr'iWi, 1^2,^3, a4,(r2,cr3), i = 2,3; (5j = (5j(/3i, /32, /33, 04, 0-2, 0-3), j = l,4. 



10 



GUI-QIANG CHEN VAIBHAV KUKREJA 



where we have omitted Ub for simphcity. 
Prom |10] . we know that 

1 



Differentiate ()2.27p with respect to /3i and let 0:4 = /3i = /32 = /^s = 0, (T2 = CT20, and 1T3 = a^Q. 
We obtain 

ri([/+) = di3^64r4,iU+) +di3-y^Gas{(73o,cr2o;U-) 

+ dj3^a'2 Ga2(f^30,'72o; f/-) V(7G(o-3o,(T2o;f^-) • ri(C/-). 

By Lemma 2.1, we have 

det (r4 ( ^7+ ) , G^^ (0-30 , 0-20 ; C^- ), (^^30 , cr20 ; f^- ), ri ( ) ) 



I "9/31 "^i I 



det(r4(C/+), G^3 (0-30, 0-20; f/-),G<x2 (0-30, 0-20; C^-),Vc/G(cr3o, 0-20; f/-) • ri(;7_)) 



K4(C/+)Ki(C/_)p?_?x?_e'^2o+'^3o(A4(C/+)e2'^2o+^3o _ \^{U-)) 
2Ki(?7+)A4(?7+)e'^2o+<T3o 



/ci(C/_)(A4(C/+)e2'^2o+<^3o + A4(?7_)) 
However, for the reflection coefficient |-K^24|j we have 



> 0. 



|5/3i(54 



det(ri([/+), G^3((T3o, 0-20; C/_ ) , G<^2 (^"so, 0-20; U^), V(/G(cr3o, 0-20; C/-) • ri(C/_)) 



det (r4 ([/+), G^3((T3o, 0-20; [/_ ) , G^j (0-30 , 0-20; f/_), V[/G(fT3o, (T20; f/- 
Ki([/„)Ki(C/+)p_u_e''2o+^3o( _ p_u_Ai(C/_) + e'"20p+u+Ai([/+)) 



K4(C/+)Ki(J7_)/9?.n^e'^2o+<^3o (A4(C/+)e2<^2o+<T3o _ Ai(C/_)) 



-A4(t/4 



g2(T2o+a"3o 



+ A4(?7_) 



< 1, 



where \K- 



24 



A4(C/+)e2'^2o+<^3o + A4(?7_) 
is not necessarily bounded away from zero, but is less than one. 



ri(f/-)) 



3. The Wave-Front Tracking Algorithm and Global Existence of Weak Solutions 

We start off here with a brief description of the wave-front tracking method to be employed 
throughout in Sections 4-7 and then establish the existence of entropy solutions when the pertur- 
bation of the incoming flow has small total variation at x = 0. 

The main scheme in the method of wave-front tracking is to construct approximate solutions 
within a class of piecewise constant functions. At the beginning, we approximate the initial data 
by a piecewise constant function. Then we solve the resulting Riemann problems exactly with the 
exception of the rarefaction waves, which are replaced by rarefaction fans with many small wave- 
fronts of equal strengths. The outgoing fronts are continued up to the first time when two waves 
collide and a new Riemann problem is solved. In this process, one has to modify the algorithm 
and introduce a simplified Riemann solver in order to keep the number of wave-fronts finite for all 
X > in the flow direction. We refer the reader to Bressan [3l [5] and Baiti-Jenssen [1] for related 
references. 



3.1. The Riemann Solvers. As seen in Section 2, the solution to the Riemann problem 
is a self-similar solution given by at most five states separated by shocks, vortex sheets, entropy 
waves, or rarefaction waves. To connect the state Ua to Ub, there exist G^-curves 77 — )• ip{rj){U) 
with arc length parametrization such that 

Ub = ^mUa) := T4(r/4) o • • • o Ti(r?i)(C/„) 

for some r? = (771, . . . ,7/4), and Uj = Tj{7]j) o • • • o Ti(r?i)(C/a), j = 1 . . . 3. 
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Next, we discuss the construction of front tracking approximations for our initial-boundary value 
problem. Let denote the initial approximation parameter. For given initial data U and with 
I? > 0, consider U a sequence of piecewise constant functions approximating U in the L^-norm, 
and the wall boundary is also approximated as described in Section 2. Set to be the total 
number of jumps in the initial data U and the tangential angle function of the wall boundary. Let 
(5^ > be a parameter so that a rarafaction wave is replaced by a step function whose "steps" are 
no further apart than 6^. The discontinuity between two steps is set to propagate with a speed 
equal to the Rankine-Hugoniot speed of the jump connecting the states corresponding to the two 
steps. At any time, the simplified Riemann solver (defined below) is employed, the constant A 
denotes the speed of the generated non-physical wave, which is strictly greater than all the wave 
speeds of system (2.1). Note that the strength of the non-physical wave is the error generated when 
the simplified Riemann solver is applied. 

Accurate Riemann solver. The accurate Riemann solver (ARS) is the exact solution to the 
Riemann problem, with the condition that every rarefaction wave {w, Rd{'w){a)}, d = 1,4, is divided 
into equal parts and replaced by a piecewise constant rarefaction fan of several new wave-fronts of 
equal strength. 

Simplified Riemann solver. When only weak waves are involved, the simplified Riemann 
solver (SRS) here is the same as the one described in [Tl[5]. That is, all new waves are put together 
in a single non-physical front, travelling faster than all characteristic speeds. In the case of a weak 
wave interacting with the strong vortex sheet /entropy wave, the purpose of (SRS) is to ignore the 
strength of the weak wave, while preserving the strength of the strong vortex sheet/entropy wave, 
and to place the error in the non-physical wave in the following manner: 

Case 1 : A weak wave {U^, Ui} collides with the strong vortex sheet/ entropy wave {Ui, [/+} from 
below. The Riemann problem U^} is solved as follows: 

for I <x(f/i,f/+), 
U2 iovx{Ui,U+)<l<\, 
U+ for I > A, 

with x{Ui, C/+) as the speed of the strong vortex sheet /entropy wave, and the state U2 is solved in 
a way that {U-,U2} is the strong vortex sheet/entropy wave starting from [/_ and x(C/i,C/+) = 
x{U-, 112)- Hence, we find that (SRS) keeps the same strength of the strong vortex sheet/entropy 
wave, and the error appears in the non-physical fronts. 

Case 2: A weak wave {U2, f7+} collides with the strong vortex sheet/entropy wave {U-, U2} from 
above. The Riemann problem {U-, [/_)_} is solved as follows: 

[/_ for I <x(f/-,C/2), 
U2 for x(C/-,f/2) < I < A, 
U+ for I > A, 

with x{U-, U2) denoting the speed of the strong vortex sheet /entropy wave. 

3.2. Construction of Wave Front Tracking Approximations 

Given the corresponding front tracking approximate solution U^{x,y) is built up as follows. 

— 19 

At X = 0, all the Riemann problems in U are solved by using the accurate Riemann solver. 
Furthermore, one can change the speed of one of the incoming fronts so that, at any time x > 0, 
there is at most one collision involving only two incoming fronts. Of course, this adjustment of 
speed can be chosen arbitrarily small. Let be a fixed small parameter with — )■ 0, as — )■ 0, 
which will be determined later. For convenience, the index j in aj will be dropped henceforward, 
and we will write Oj as a when there is no ambiguity involved; the same applies for /?; and we will 
moreover employ the same notation a as a wave and its strength as before. 

Case 1: Two weak waves with strengths a and j3 interact at some x > 0. The Riemann problem 
produced by this collision is solved in the following way: 
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• If \af3\ > u}^ and the two waves are physical, then the accurate Riemann solver is employed. 

• If \a/3\ < uj^ and the two waves are physical, or there is a non-physical wave, then the 
simplified Riemann solver is employed. 

Case 2: A weak wave a interacts with the strong vortex sheet/ entropy wave and one weak wave 
at some x > 0. The Riemann problem produced by this collision is solved in the following way: 

• If |q| > and the weak wave is physical, then the accurate Riemann solver is applied. 

• If |a| < uj-o and the weak wave is physical, or this wave is non-physical, then the simplified 
Riemann solver is applied. 

Case 3: The flow perturbation due to the Lipschitz wall boundary. 

• When the change of the angle of the boundary is larger than and the weak wave is phys- 
ical, then the accurate Riemann solver is employed to solve the lateral Riemann problem. 

• If the change of the angle of the boundary is less than u^, then this perturbation is ignored. 
Case 4-' The physical wave collides with the boundary. The accurate Riemann solver is employed 

to solve the lateral Riemann problem. 

Case 5: The non-physical wave collides with the boundary. We can allow these waves to cross 
the boundary. 

3.3. Climm's Functional and Wave Interaction Potential 

The goal in this subsection is to construct the suitable Glimm-type functional and the associ- 
ated wave interaction potential Q for our initial-boundary value problem. This involves a careful 
incorporation of the additional nonlinear waves generated from the wall boundary vertices. 

Definition 3.1 {Approaching waves), (i) Two weak fronts a and /3, located at points ya < yp 
and of the characteristic families j^, jp S {1, . . . ,4}, respectively, are said to be approaching each 
other if the following two conditions are concurrently satisfied: 

• ya and yp are both in one of the two intervals into which M is partitioned by the location 
of the strong vortex sheet /entropy wave. That is, both waves are either in VL^ or 

• Either ja > jp or else ja = jp and at least one of them is a genuinely nonlinear shock. 
In this case, we write (ce, j3) G A. 

(ii) We say that a weak wave a of the characteristic family ja is approaching the strong vortex 
sheet /entropy wave if either a G 0_ and ja = 4, or a G and ja = 1. We then write a G ^t,/e. 

(iii) We say that a weak wave a of the characteristic family ja is approaching the boundary if 
a G il_ and ja = ^- We then write a G Ab. 

Define the total (weighted) strength of weak waves in U^{x, •) as 

V{x)=^\ba\. 

a 

Here, for a weak wave a of the j-family, its weighted strength is defined as 

I k^a if a G ^2+ and ja = 1, 
ba = \ 

la if a G 

where k^ = "^^r and the coefficient K21 as given in Lemma 2.6. 
Next, the wave interaction potential Q{x) is defined as 

Q{x) = C* \bahp\+K* K\+Y\^P\+^Y.\'^^ 

= QA + Qv/e + Qb + Qe- 

Here the constants K* G (i^ii, 1) and A'^q > K^o^ while C* is a constant to be specified later. To 
control the total variation of the new waves produced by the boundary vertices, Qe in our wave 
interaction potential Q{x) is an added term, compared to that for the Cauchy problem. 



(3.1) 



(3.2) 
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The Glimm-type functional Q is defined as follows 

g{x) = V(x) + KQix) + \U^{x) -U+\ + |C/o(x) - C/o"|, (3.3) 

where the states U<^{x) and C/^(x) are the below state and the above state of the strong vortex 
sheet/ entropy wave respectively at "time" x, Uq and Uq are the below and above state of the strong 
vortex sheet/entropy wave respectively at x = 0, and k is a large positive constant to be determined 
later. 

We remark that V, Q, and Q remain unchanged between any pair of subsequent interaction times. 
However, we will demonstrate that, across an interaction "time" x, both Q and G decrease. 

Proposition 3.1. Assume that TV(C/o(-)) + TV(5' (•)) is sufficiently small. Then V{x) will remain 
sufficiently small for all x > 0. Moreover, the quantity TY{U^ (^j ")) ^'^■^ cl uniform bound for any 
t9 > 0. 

Proof. With the Glimm-type functional Q, consider 

Ag{x) = g{x+)-g{x-), 

where x~ and x"*" denote the "times" before and after the interaction "time" x > 0, respectively. 

Case 1: Two weak waves a and f3 collide. Then the states (x) and Uo (x) do not alter across 
this interaction "time" x > 0. Hence, we have 

Ag{x) = V(x+)-V(x-) + k(Q(x+)- Q(x-)) 

< Bi\bab^\ + K {-C*\babp\ + C*Kb^\Vix-) + Bo\bab^\) , 

where Bq and Bi are constants independent of -d. 

Case 2: A weak wave a of the 1-family interacts with the boundary. 

Ag{x) = Kbia -a + K {C* Ki,iV{x^)a + K*Khia - a) , 

where K*Ki,i < 1. 

Case 3: A new 4-wave a produced by the Lipschitz wall boundary. 

Ag{x) = Kboei + K (c*Kbo0iVix-) + K*Kboei - K^oOi) , 

where Ki^q < K^o is large. 

In the following two cases, the states Uoix) and U^{x) change across this interaction "time" 
X > 0. 

Case 4^ A weak wave a of the ^'fo-'^i^y collides with the strong vortex sheet/ entropy wave from 
below. 

Ag{x) = V(x+) - V(x~) + |[/^(x+) - [/^(x-)| 

+|C/o(x+) - f/o(x-)| + K (Q(x+) - Q(x-)) 

4 

= ^^Kija - a + K(^C*{KiiV{x~)a + Ki4V{x~)a) - K*a + KiiQ^. 
i=i 

Case 5: A weak wave a of the 1-family collides with the strong vortex sheet/entropy wave from 
above. 

4 

Ag{x) = K2ja -ba + k(c* {K2iVix-)a + K2AVix-)a) - K*b^ + i^ai") • 
i=i 

In the cases above, Kn < K* < 1, b^ > 2X2x101 in connection with the weight k^, and the 
constant C* > ^Bq > is large. 

Next, we establish that the total (weighted) strength of waves in U'^{x, •) remain sufficiently small 
for all X > if it is sufficiently small at x = 0. More precisely, 

V(x) <C 1 for ah x > 0. 
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This can be proved as follows: 

(i) "Time" xi > is the first interaction. Given that V (xj") = V(0) < TV(t^(-)) < 1 and 
Yl'iZo ^ TV(g''(-)) ^ 1 in Cases 1-5 above, we conclude that, for k sufficiently large and lo^ small 
enough, 

Ag(xi)<o, i.e., g{xi) <g{x^) = g{o). 

Therefore, 

v(x+) < Gixi) < g{o) < v(o) + kQ{o) 

oo 
1=0 

oo 
1=0 

(ii) V(x~) <C 1 and ^(x^) < ^(x~) /or any m < n. Then, for the next interaction "time" Xn, 
similar to Case 1, we also conclude 

Ag(x„)<o, i.e., g{xt)<g{x~) = g{x+_,). 

Therefore, all together, we obtain 

V(x+) + |[/o(x+) - C7+I + \U,{xt) - 

<g{xt)<G{x-) = g{xt_,)<...<g{o) 

= V(0) + kQ(0) 

oo 

= V(0) + Hi(c*V\0) + V(0) +KboY.^i) 

1=0 

oo 
1=0 

This implies that V{x) <^ 1 for all x > 0, since C is independent of x. 

Furthermore, the total variation of U^{x,-) is uniformly bounded. More precisely, we conclude 
that 

TV{C/''(x, •)} « V{x) + |?7^(x) - f/o+l + |C/o(x) - ^7o-| + |fT2o| + kaol = (3.4) 
This completes the proof. 

In order to have a front tracking approximate solution U^{x, •) defined for any time x > 0, along 
with a uniform bound on the total variation, we also need to have that the number of wave-fronts 
in U^{x, •) is finite. This is given by the subsequent lemma. 

Lemma 3.2. For any fixed "i? > small enough, the number of wave- fronts in (x, y) is finite 
and the approximate solutions (x,?/) are defined for all x > 0. Moreover, for any x > 0, the 
total strength of the all non-physical waves is of order 0{1) (5^ + w^). 

Proof. We first note the total interaction potential Q(x) remains unchanged when there is no 
interaction and decreases across an interaction "time" x > as discussed in Cases 1-5 in Proposition 
3.1. Furthermore, from Cases 1-5 and the subsequent analysis above, we have concluded that 
V(x) <^ 1. Hence, one can fix some number v £ (0, 1) such that 

AQ(x) = Q(x+)-Q(x-) 

{—I'lbabpl if both waves a and /3 are weak, 
— i^|5q,| if the weak wave a hits the strong vortex sheet /entropy wave, (3.5) 
— i^l^il if the angle of the boundary changes. 

Now, following an argument similar to the one given in [1], we reach the following conclusions. 
Note that initially Q(0) is bounded and Q decreases thereafter for each case. Moreover, in the 
case where the interaction potential between the incoming waves or the change of the angle of the 
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boundary is larger than w^, Q decreases by at least voj^ in these interactions, as implied by the 
bounds given in (|3.5p . Following the wave-front tracking method in our problem, new physical 
waves can be only produced by such interactions. Furthermore, when the weak wave a of 1-family 
collides with the wall boundary, we solve the lateral Riemann problem and have shown earlier that, 
after this interaction, there is only a reflected wave of 4-family with the reflection coefficient 1. 
Hence, before and after this interaction, the number of the waves keeps the same, and this implies 
that the number of the waves is finite. Finally, because non-physical waves are generated only 
when physical waves collide, we can also conclude that the number of non-physical wave fronts are 
finite; and, provided that two waves can only collide once, the number of interactions is also finite. 
Consequently, it follows that the approximate solutions U'^{x,-) are defined for all times x > 0. 
The similar argument allows us to conclude that the total strength of all non-physical wave fronts 
at any x is of order 0{1){5^ + ui^). This completes the proof Lemma 3.2. 

Following the line of arguments given in [H [3] for the wave-front tracking algorithm and Lemma 
3.1 above, we finish this section with the following theorem for the global existence of entropy 
solutions to the initial-boundary value problem (jl.ip and ()1.7p - ()1.8p . 

Theorem 3.1. Suppose that TY {Uo{-)) -\-TV {g' {■)) is small enough. Then, for the initial-boundary 
value problem (jl.ip and ()1.7p - (jl.Sp . there exists a global weak solution in BV satisfying the steady 
Clausius entropy inequality ()1.9p . 

4. The Lyapunov Functional for the -Distance between Two Solutions 

To show that the front tracking approximations, constructed for the existence analysis in Section 
3, converge to a unique limit, we estimate the distance between any two ^^-approximate U and V of 
problem (jl.ip and ()1.7p - ()1.8p of initial-boundary value type. To this end, we develop the Lyapunov 
functional (^{U,V), equivalent to the L^-distance: 

\\u{x, ■) - v{x, oiili < y)<c \\u{x, ■) - v{x, oiiii, 

and prove that the functional <^(C/, V) is almost decreasing along pairs of solutions, 

^ {U{x2, •), V{x2, ■)) - ^ {U{xi, ■),V{xi, ■)) < C^{x2 - xi), for all X2>xi> 0, 

for some constant C > 0. Here U and V are two approximate solutions constructed via the wave- 
front tracking method, and the small approximation parameter i? is responsible for controlling the 
subsequent errors: 

• Errors in the approximation of the initial data and the boundary. 

• Errors in the speeds of shock, vortex sheet, entropy wave, and rarefaction fronts. 

• The total strength of all non-physical fronts. 

• The maximum strength of rarefaction fronts. 

Along the line of arguments presented in [HI [211 [23], with "time" x fixed, at each y, one connects 
the state U{y) with V{y) in the state space by going along the Hugoniot curves ^i, C2, C3, and 5*4. 
Depending on the location of the strong vortex sheet/entropy wave in U (y) and V{y), the distance 
between U{y) and V{y) is estimated along discontinuity waves in possibly different "directions", 
determining the strength of the j-Hugoniot wave hj{y) in the following way: 

• Suppose that U{y) and V{y) are both in r2_ and Then one begins at the state U{y) 
and moves along the Hugoniot curves to reach the state V{y)- 

• Suppose that U{y) is in and V{y) is in Then one begins at the state U{y) and 
moves along the Hugoniot curves to reach the state V{y). 

• Suppose that V{y) is in r2„ and U{y) is in Then one begins at the state V{y) and 
moves along the Hugoniot curves to reach the state U{y). 
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Define the L^-weighted strengths of the waves in the solution of the Riemann problem (U (y), V{y)) 
or (y{y),U{y)) as follows: 



hj{y) whenever U{y) and V{y) are both in Vl^ 



qj(y) = < u;™ • hj{y) whenever U{y) and V{y) are both in different domains, (4.1) 
Wj ■ hj{y) whenever U{y) and V{y) are both in 0+, 

with the constants w'j, wj^, and Wj above to be specified later on, based on the estimates of wave 
interactions and reflections in Lemmas 2.2-2.7. 
We define the following Lyapunov functional, 

4 oo 

^{U,V) = Y, J hivWMdy, (4.2) 

where the weights are given by 



Wj{y) = 1 + K^Ajiy) + K2 {Q{U) + Q{V)) . (4.3) 

The constants ki and K2 are to be determined later. Here Q denotes the total wave interaction 
potential incorporating the boundary effect as defined in ()3.2p . and Aj{y) denotes the total strength 
of waves in U and V, which approach the j-wave qj{y), defined in the following manner (for y where 
there is no jump in U or V): 

A / \ ^ / N ^ , \ { Hjiy) if 7-wave qAy) is small and the v'-field is genuinely nonlinear, 
AJy) = FAy) + GAy) + { ^ l!^^^' / n „■ , 

3\yj 3\yj < -^^^^ I if J = 2, 3 and qj{y) = B is large. 

(4.4) 

Next, we define the following global weights Gj: 



GM = 


U, V are both in 


U, V are in distinct regions 


U, V are both in 


Gi{y) 


4B 


2B 


4B 


G2,M 











GM 


4B 


2B 


2B 



Under the assumption that TV(C/o(-)) + TV(Vo(-)) + TV(5'(-)) is small enough with U{x, •), V{x, •) 
G BV n L^, one concludes 



M-'\\U{x,-)-V{x,-)h^<J2 



h 



\qj{y)\dy<M\\U{x,-)-V{x, 



l<Wjiy)<M, j 

where the constant M is independent of and "time" x. Here we define the strength of any large 
wave of the 2- or 3- characteristic family to equal to some fixed number B (bigger than all strengths 
of small waves), and the concepts "small" and "large" mean the waves that connect the states in 
the same or in the distinct domains and respectively. 
The summands in (14.41) are defined as follows. 




E 



a&J\SC / 
ya>y,i<ka<j 



iJ2a&J{U)\SC,yc<y,k^=j^^a€j{V)\SC,yc>y,kc=j)\^\ 9j(y) < 0, 

(J^aCiJ{V)\SC,ya,<y,k^=j + ^aCiJ{U)\SC,y^>y,k^=j)\'^\ ^jiv) > 0' 

where, at each x, a stands for the (non-weighted) strength of the wave a G J', located at the point 
ya and belonging to the characteristic family ka] J = J{U) U JiV), SC = SC{U) U SCiy) is the 
set of all waves (in U and V) and the set of all large (strong) characteristic discontinuities (in U 
and V) respectively. 
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Consequently, there holds 

C~^\\U{x, •) - V{x, OIIli < V) < C\\U{x, •) - y(x, OIIli , (4.5) 

for any x > with the constant C > depending only on the quantities independent of x: the 
strength of the strong vortex sheet/entropy wave and TV(C/o(-)) + TV(Vb(")) + TV(5'(-)). 

We now analyze the evolution of the Lyapunov functional in the flow direction x > 0. For j = 
1, . . . , 4, we call \j{y) the speed of the j-wave qj{y) (along the Hugoniot curve in the phase space). 
Then, at a "time" x > which is not the interaction time of the waves either in U {x) = U (x, •) or 
V{x) = V{x, •), an explicit computation gives 

±<p{U{x),V{x)) 
ax 

4 4 
4 

= E E )|WG(ya ) {vo. - A,(y^)) - \qMWAyt) - Hvt))) 

a&J j=l 
4 

+ ^\qMWj{b){m + Xj{b)), (4.6) 
i=i 

where i/a denotes the speed of the Hugoniot wave a £ J', b = g{x)'^ stands for the points close to 
the boundary, and yb is the slope of the boundary. 
We present the notation 



Ea, = \ql\W; -y^)- \g-\Wr ^XJ - y„j , (4.7) 
Eb,, = \qj{b)\W,{b){yb + \j{b)), (4.8) 



where qf = q,{yt), Wf = W,{yt), and Xf = A,(y±). 
Then (14.611 can be written as 



d 4 4 

-^{U{xlV{x)) = Y.Y.E-,+Y.E^,r (4-9) 

aGj j=l j=l 

Our central aim is to prove the bounds: 

4 

^Eo,j <0{l)'&\a\ when a is a weak wave in J , (4.10) 
i=i 

4 

^Eaj<0{l)\a\ when a is a non-physical wave in J , (4.11) 

4 

Eoij < when a is a strong vortex sheet/entropy wave in J ^ (4.12) 

4 

E'fcj- < near the boundary, (4.13) 

where the quantities denoted by the Landau symbol 0{\) are independent of the constants ki and 

From ()4.10p - ()4.13|) together with the uniform bound on the total strengths of waves ()3.4p . we 
obtain 

^^iU{x),V{x)) <0(l)^. (4.14) 
ax 
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Integration of (|4.14p over the interval [0, x] yields 

^ {U{x), V{x)) < <P (C/(0), y(0)) + 0{l)Sx. (4.15) 

We remark that, at each interaction "time" x when two fronts of U or two fronts of V interact, 
by the Glimm interaction estimates, all the weight functions Wj{y) decrease, if the constant K2 
in the Lyapunov functional is taken to be sufficiently large. Furthermore, due to the self-similar 
property of the Riemann solutions, decreases at this "time". 

In the next section, we establish the bounds (|4.1U|) - (j4.13|) . particularly (|4.1'2|) and (|4.13|) . when 
a is a strong vortex sheet /entropy wave in J and near the Lipschitz wall boundary, respectively. 



5. The L^-Stability Estimates 

For the case of the non-physical waves in ^7, as well as the case that the weak wave q G ^7 := 
J{U) U J{V), which appears when U and V are both in or estimates ()4.10p and ()4.1ip 
are shown similarly based on the arguments in Bressan-Liu-Yang 0, provided that |^^^^^_^^^^| is 
sufficiently small and ki is sufficiently large. In what follows, we focus only on the other two cases, 
namely ()?l2D and (1^3]) . 

Case 1: The first strong vortex sheet/ entropy wave a in U or V is crossed. Then, by Lemma 
2.4, we have the estimates: 

/i+ = h-^+Kiihl, (5.1) 

ht = Kuh^. (5.2) 

Moreover, the essential estimate < 1 given in Lemma 2.4 ensures the existence of desired 

weights w\ and in the following way. 

Lemma 5.1. There exist w\, w\, and jt, satisfying 



w 



4 < 1, (5.3) 



w 



—All 



X 



1 — A2,3 



-^4 — -^2,3 



< 76 < 1- (5.4) 



With Lemma 5.1, we estimate Ej for j = 1, . . . ,4, starting with Ei: By (j5.ip and ([57 

El = \q^\{\^ -ijaXW^ -Wr) + W^{\qt\{\t -ya)-\qi\{Xi -Va)) 
= 2Ski^|/i^||A5~ - yd + {\qT\\K - Val - \qt\\^i - 

< 2BKiw^^\h^\\X^ - ijal + W+ (w\\ht\\X^ - y„| + w'lKu\h^\\X^ - y„| - wT\ht\\Xt - i/a 

< 2BKiw\\h^\\X^ - ya\ + 2Bki (^w\\hf\\X'^ - ya\ + 7f)^^4l^4 1(^4 - Va) - w"'-\hf\\Xf - y„ 
+ {kiA^+ + M)\q^\\X^ - ijal - {KiA^++M)\qt\\Xt -ya\, 

where = Wi{y+) = 2Bki + HiiAy^+ + M , and M = 1 + K2{Q{U) + Q{V)) is a positive constant. 
The term A-^^+ = Fi(y+) + Hi{y^) here is the total strength of all the weak waves in U and V 
which approach the 1-wave qf = qi{ya)i a-^d the term 2Bki is from the weight Gi(y^). 
For j = 2, 3, since Wj~ = WJ' , (|4.6p reduces to 

Ej=W-{\q'^\{X+ -y^)-\q-\{X- -y^)) <0{1)b[^+ E + 

j^{2,3} 

where k 7^ {j, 1,4}. 
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For j = 4, 

Ea = Iq^liK -ya)iWl -W^) + Wl{\qt\{Xt -ya)-\qA\{K -Va)) 
= -2BKi\q^\{X^ - yo.) + {2Bki + KiA^+ + M)\ql\{Xl - i/a,) 

4 

-{2Bki + KiA^+ + M)\q^\{X^ - y^), 

where = Wiiy^) = 2Bki + KiAy^+ + M , and M = I + K2{Q{U) + Q{V)) is a positive constant. 
The term Ay^+ = Fi{y'^) + Hi{y^) is the total strength of all the weak waves in U and V which 

approach the 4- wave q'l = ^4 (?/+), and the term 2Bki is from the weight G4(y+). 

For the weighted L^-strength qj{y) in (|4.ip . we choose w\ small enough relatively to w"'' and w\ 
large enough relatively to w'f'^ choose ki large enough and the total variation of U and V so small, 
and use (I5.1D-(I5.2I) to obtain 



4 

Y,Ej < 2BKi(w\\ht\\X^ -y^\+-fbwl\h^\iX^ -ya)-wT\ht\\Xt -ya^ 



+{'iiA^+ + M)\qY\\X^ - y^\ - {kiA^+ + M)|g+||A+ - 

+2BKiw'i\h^\\X^ -ijal- 2Bni\q^\{X^ - y„) 
+ {2Bki + KiA^+ + M)w'l'\Kup^\{Xt - y^) 

-{2Bki + KiA^+ + M)\q^\{X^ - ija) 
+0(1)5(2^ + 2 \q7\ + {\Q2\ + \13 



i^{2,3} 

-{1 - -fb)2BKiw'l\h^\{X^ - ija) 

+{f^iA^+ + M){wT\KuK\{Xt - i/o,) - wl\h^\{X^ - i/a,)) 

+2BKiw'li\ht\ + |/irl)|Ar - ya\ - 2BKiwT\ht\\Xt - y«| 

+2BKiwf\KuK\{X+ - i/a) - 2BKiwl\h^\{X^ - y^) 

+(/ti^^+ + M)w\\hY\\X^ - ya\ - {kiA^+ + M)t/;f |/i+||A+ - i/a 

+0{1)b(2^ + 2 \q7\ + {\<12\ + \% 



i^{2,3} 

< 0. 

Case 2: The weak wave a between the two strong vortex sheets /entropy waves in U and V is 
crossed. For j = 1, we have 

El = \qt\{W+-Wr){Xf-y^) + W^{\qt\{X+-y^)-\qT\{X^-y^)) 

< Hii\qf\\a\\Xf - ya\ + 2BHii ((|g+| - \ql\){Xt - ija) + |grl(A+ - A^)) 

< K^\qt\\a\\Xf - y„| + 2Bk^ ((|g+| - \q'^\){Xt - ija) + 0{l)\q^\\a\) . 

For the cases when j = 2, 3, we have 

E, = B[{w;-Wr){xf-y^) + W^{Xf-Xj) 
< B(^-Ki\a\\X+ -y^\+0{l)\a\ 

For j = 4, we have 

Ea = \qt\iWl-W^){Xf-ya) + W^{\qt\{X+-y^)-\q^\{X^-y^)) 

< '^ikj l|a||Af - ya\ + 2Bki - \q^\)iX+ - ij^) + \q^\{Xt - Aj)) 

< ^^l\qt\W\\>i - + 2i?A^i {{\qt\ - \q^\){Xi - ijo) + 0{l)\q'^\\a\) . 
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Then we have 

4 

^E,<^,0{l)[-\a\ + \a\ J2 \^k\ + \Qk\+ E k,t| - ICl) + 
Since — |g^r|| < — < 0(l)|a| when k ^ {2,3}, we obtain 

4 

if all the weights vf^ are small enough and k\ is sufficiently large. 

Notice that the choice of the upper or lower superscripts depends on the family number k^. 

Case 3: The second strong vortex sheet/entropy wave a in U or V is crossed. For this case, by 
Lemma 2.6, we have 

= K2iht, (5.5) 
/i4 = h+ + K24ht. (5.6) 

Moreover, the essential estimate |i^24| < 1 in Lemma 2.6 ensures the existence of desired weights 
vuf and vul in the following manner. 



Lemma 5.2. There exist wf, w1, and ja satisfying 



w 



-^4 — -^2,3 



i^24 < 7a < 1. (5.7) 



^1 — A2,3 

With Lemma 5.2, we estimate Ej for j = 1, ... ,4 as follows: By (|5.5p . 

El = \q^\{X^ -ya){W^ -W^)+W+{\qt\{Xt -ya)-\qi\{Xi -i/a)) 
= -2BKi|gf ||A^ - ijal + (4Ski + kiA^+ + M)\q+\{X+ - y„) 

+ wT\K2iht\{4BK^ + + M)\XY -ijal, 

= -2BKi\qY\\X^ -y^\+wT\K2ih+\{4BKi + KiA^++M)\X^ -y^\ 

-uj1\h+\{2BKi + KiA^+ + M)|A+ - ijal - 25ki<|/i+||A+ - 

where = Wi{y+) = ABki + KiAy^+ + M , and M = 1 + K2{Q{U) + Q{V)) is a positive constant. 
The term Ay^,+ = Fi[y^) + Hi{y^) here is the total strength of all the weak waves in U and V 
which approach the 1-wave qf = (/i (?/„), and the term ABki is from the weight Gi{y^). 
For j = 2, 3, since Wj~ = , (|4.6|) reduces to 

Ej = W+{\q+\{X+-y^)-\q-\{X--y^)) 

< 0{1)b{{}+ ^ \qt\ + \qt\): 

where k ^ {j, 1,4}. 
By §M) and 

Ea = Iq^m -ya){W+ -W^) + W+{\q+\{Xt -ya)-\qI\{X^ -Va)) 

= {\qt\{Xt - ya) - krK^r - Va)) 

< W+ {w2\h^\iXi -y^) + w2K24\h+\{Xi -ya^)-wT\K\{X^ -ya)) 

< 2Bki {w2\h^\iXt - ija) + laW'i\h't\\X+ " l/a | " < | /l^ I (Aj " Va)) 

+ {KiAy^+ + M)\qt\{Xt - ija) - (kiA + +M)|g4 KAJ -y«), 
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where = Wi{y+) = 2kiB + kiA^+ + M , and M = 1 + K2{Q{U) + Q{V)) is a positive constant. 

The term = F^d/^) + Hi{y'^) here is the total strength of all the weak waves in U and V 

which approach the 4- wave = (74(7/^), and the term 2Bki is from the weight G4^{y'^). 

For the weighted L^-strength qj{y) in (4.1), when w1 is small enough relatively to w'^ and wl 
is large enough relatively to w'f' , ki is large enough, applying (|5.5p and ()5.6p . the total variation of 
u and V is so small that 
4 

Y,Ej < 2SKi(w;2|/irl(A|-ya)+7a<|/i|||A^-?/„|-<I^II(Ar-ya)) 
3=1 

+ {kiA^+ + M)\qt\{\t - ija) - {t^iA^+ + M)\q-^\{Xl - ij^) 

-2BKi\q{\\X{ - ya\ - 2BKi\q^\\Xl - ya\ 
+ {ABki + KiAy^+ + M)wX'\K2ih'[\\\^ - y^l 

-{2Bki + KiAy^+ + M)w'i\h'[\\Xl - ya\ 

+0{l)B(2d + 2 + (1^+1 + 1^3+1) 

i^{2,3} 

= -{l--fa)2BKiwt\h+\\Xt -ya\ 

+(^1^1^+ + M)[wT\K2iht\\X^ - y„| - wl\ht\\\t - ya\) 

+2BKiwl\hl\{Xl - ijc) - 2BKiw'^\hl\{\l - y^) 

+2BKi{-w'{'\K2iht\ + 2w'{'\K2ih+\)\X^ - y^\ - 2BKiwl\h+\\Xl - y^\ 
+{kiA^+ + M)wl\h+\{Xl - ija) - {kiA^+ + M)w'X\hl\{\l - ijo:) 

+0(1)5(2^ + 2 Yl \^t\ + {\^t\ + \4\)) 

i5^{2,3} 

< 0, 
which yields (|4.10p . 

Case 4: Close to the Lipschitz wall boundary. This case differs from the Cauchy problem. 
Here we will use the particular property of the boundary condition (jl.Sp : The flows of U and V 
are tangent to the Lipschitz wall, implying that they must be parallel with each other along the 
boundary. Then a piecewise constant weak solution is constructed only along the Hugoniot curves 
determined by the Riemann data U{b) and V{b), the states of solutions U and V, respectively, close 
to the boundary. 

Proposition 5.2. Suppose that U{b) = {u,v,p,p) and V{b) = {u,v,p,p) are states in a small 
neighborhood 0£{U-) of U- satisfying ^ = ^ = Zh, and v,v ^ 0. Denote by hj{b) the strength of 
the shock in the Riemann problem determined by U{b) and V{b), and by Xj the corresponding 
j^^ -characteristic speed. Then 

|Aj -ifcl ~ |/ii(6)|, i = 2,3, (5.8) 
1/14(6)1 < \hi{b)\ + 0{l)\h2mX2 - Zb\ + |/ii(6)|0(l)|ife|, (5.9) 

\h,{b)\=0{l)\h,{b)\, 1<0(1)<| (5.10) 
where z^ is the slope of the Lipschitz wall. 

Proof. We do the proof by analyzing the following two cases. 

Case 1: hi{b) = and h^{b) = that corresponds to the case p = p. Starting at the state Ub, we 
move along the Hugoniot curves of the second and third families to reach Vb- Note that these two 
families are the contact Hugoniot curves, and so A2 and A3 are constant along the Hugoniot curves. 



22 



GUI-QIANG CHEN VAIBHAV KUKREJA 



Given that A2,3 = ^, r2 = (1> 0) 0)"'', and r-^ = (0, 0, 0, p)"*", the quantity ^ remains unchanged 
as the initial value , i.e., Zf, in this process by the boundary condition (jl.Sp . Therefore, we 
conclude that A2,3 = ifc, equivalently, 

4 - A2,3 = 0. 

Case 2: hi{b) ^ that corresponds to p ^ p. Starting at the state U{b), we move along the 
1-Hugoniot curve to reach Ui, then possibly move along the 2-contact Hugoniot curve to reach U2, 
the 3-Hugoniot curve to reach U3, and the 4-Hugoniot curve to reach V{b). 

To make clear some essential relations among the strengths /i2(&), /i3(fe), and /i4(6), we 

project {u,v,p,p) onto the {u,v)"pla,ne. Let rij^ be the projection of ri onto the u-axis, r2\(^u,v) be 
the projection of r2 onto the (u, f)-plane; and so on. At the background state C/_, there holds 

^i\u — ~T^i\ui ri|j, = r4|t,, ril(p^p') = ~r4|(p,p)) 

r2 = r2 1 TC3\{u,v) = 0- 

We first note that /i4(6) 7^ 0. Given that ri|(„^,) = ki{—Xi,l)~^ along with finite characteristic 
speeds Ai and Zb ~ 0, there always holds if, < — ^ near the state U-. So we can conclude that, 
in the (n, u)~plane, the derivative ^ along the 1-curve is always larger than Zf,- This implies that 
E(ffA^E(ff>^ Meanwhile, we have ^ = ^ = ^ and ^^4 = 1^- Hence, 

u{Ui) I u{Ui,) ' u{Ui) u{U2) u{U3) u{Vt) u(Ut) ' 

v{Ui) _ V{U2) _ V{U3) v{Vb) 



u{Ui) u{U2) u{U^) ' u{Vb)' 
Thus, we conclude that there is some distance along the 4-Hugoniot curve to reach V^- Therefore, 
/14 / 0. 

Next, we present an essential estimate to bound \ht\ more precisely in terms of |/i4|. To that 
end, define the signed length of {JJ\ — Vh)\(u,v) and (V;, — C/3)|(m,i,) by d\ and on the (u, u)-plane: 

r ii(f/i-t/6)i(„,.)ii if/ii>o, 

\ -Wi-'Ub\u^,^\\ if /ii <o. 



C?1 



and 



^4 



II(V6-c/3)Ik.)II if/i4>o, 

- f^3)|(n,,;)ll if /14 < 0. 



Secondly, we note that 



\\2-Zb\=0{X)\dx\=0[\)\hx{h)\. 

Moreover, since A2 = = = A3, we can similarly conclude 

|A3-ife| = 0(l)|di| =0(l)|/ii(6)|. 
Using the following projections on the (u, f )-plane, 

ri|n = ~r4|^t, ril^, = r4|t,, 

r2 = r2|(„,t,), r3|(„^„) = 0. 

Thirdly, we note that 

-d4 = 0(l)/i2(6)(A2-i6) + d, 
where dcost^i = diC0S(^2) V'l denotes the angle between (l,ife) and r4|(„^^), Lp2 denotes the angle 
between ri|(„ and (1, i^), = 992 + 2a for a = arctan(i6), and 

J C0SLP2 , cos((y9i — 2a) cos (^1 cos(2q) + sin 991 sin(2a) 
d = d\ = d\ = 01 - 



cos Lpi cos Lpi cos ipi 

= di (cos(2a) + C'(l)sin(2a)) =4(1 + C'(l)a) = di{l + 0(1)4) . 

Hence, there holds 

-d4 = 0{l)h2{b){X2 - Zb) + di{l + 0{l)zb). 
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At [/_, there also holds ri|(„pp) = -r^l ^u,p,p) ^-iid ^i\v = T^ilv, which implies 

di di 
hi /14 

Hence we note the following key estimate: 

- h^b) = 0(l)/i2(6)(A2 - 4) + hi{b){l + (5.11) 
Estimate (jS.lip now implies 

\h4b)\ < \hi{b)\+0{l)\h2mX2-Zb\ + \hi{b)\0{l)\ib\ 

< \hi{b)\+0{l){\h2{b)\ + \zt\)\hi{b)l 



yielding 



\hiib)\ = d{i)\h^ib)\ withi<o(i)<| 



given that |/i2(6)| + |ife| is always small enough, and this is guaranteed by the sufficiently small total 
variation of the initial perturbation Uq and the perturbation of the boundary. □ 

We note that the requirement ^ = | = i;, in Proposition 5.2 is just the boundary condition (|1.8p 



because Zj, here is the slope of the Lipschtiz wall. 

Applying Proposition 5.2 now yields 

^6,1 = \qi{b)\Wi{b){it + Xi) 

= -4BKiw\\hi{b)\\Xi\+0{l)\hi{b)\ 
= -4BKiw\\hi{b)\\Xi\+0{l)\h^{b)l 

Ebj = WMWjmzb + A,) = \qj{b)\Wj{b){ib - Xj) + 2X,\qj{b)\Wj{b) 

= 0{l)w';\hj{b)\{zb-Xj) + 0{l)Xjw''j\hj{b)\ (5.12) 
= 0{l)hi{b) = 0{l)h4{b), j = 2,3, 

Eb,4 = \qi{b)\Wi{b){zb + X^) 

= ABKiwi\h4ib)\\Xi\ + Oil)\h4{b)\ 

< ABKi\Xi\w'l{\hi{b)\ + 0(l)|/i2(6)||A2 - zb\ + 0{l)\hi{b)\\zb\) + 0{l)\h^{b)\. 
Using Lemma 5.1, we can choose Wi and w'l such that 

wl < wl 

Then, with the total variation of the incoming flow perturbation and the boundary perturbation 
small enough and ni large enough, one has 

4 

Y,Eb,j = ABKi{w'^-w\)\hi{b)\\Xi\ 
i=i 

+0{l)BKi\Xi\hl{\h2{b)\ + \ib\)\hi{b)\+0{l)\h^{b)\ 
< 6{l)4BKi{wi-w'l)\h4{b)\\Xi\ 

+0{l)BKi\Xi\wl{\h2{b)\ + \ib\)\hi{b)\ + 0(l)|/i4(6)| < 0, 

provided that |/i2(^)| + \zb\ is sufficiently small. This is guaranteed since the total variation of the 
incoming flow perturbation and the boundary perturbation are sufficiently small. 
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6. Existence of a Semigroup of Solutions 

As a corollary of the essential estimates in Sections 3-5, we can now establish the existence of 
the semigroup S generated by the wave-front tracking method, as well as the Lipschitz continuity 
of S. 

Proposition 6.1. Suppose that TV{Uq{-)) + TV{g'{-)) is small enough. Then the map 

produced by the wave-front tracking algorithm is a uniformly Lipschitz continuous semigroup satis- 
fying the properties: 

(i) §tU_^ U, _St,^t;U =_§t,+x,U, for ah xi,X2 > 0; 

(ii) \\§iU - S^FIIli <C\\U- V\\li + C^x, for ah x > 0. 

Proof. Since S'' is generated by the wave-front tracking algorithm, property (i) is immediate. 
Next, property (ii) is proved as follows. Take a pair of front tracking -(^-approximate solutions U"^ 
and of (HIT]) and (fTTll - P^ with U{-) and F(-) as the initial data, respectively. Using 
and (|4.15|) . at any x > 0, we have 

\\U^{x)-V^{x)\\li < C<^{U^{x),V^{x)) 

< C<!>{U^{0),V^{0)) + Cux 

< C\\U -V\\li +C^x. 

Hence, the ?9-semigroup is Lipschitz continuous. 
Definition 6.1. For a given i^q > 0, we define the domain: 

the set consisting of points U : R i— ?• M'* 
such that there exists one point G M and 
V = cH , fc/_, g{x)<y<y\ 
U+, y' < y, 
- U G Li(M;M^) and TV(C/ - U) < i^q. 

Remark 6.1. Given a solution U{x,y) to the initial-boundary value problem of (jl.ip and ()1.7p - 
(|1.8p . we note that, if U^{y) := U{x,y) £ V at any fixed x > 0, then > ^(O) = at x = and 
> g{x) for x > as a strong vortex sheet /entropy wave is present. 

The semigroup § generated by the wave-front tracking algorithm is provided by the next theorem. 

Theorem 6.1. Suppose that TV{Uq{-)) -\-TV{g'{-)) is small enough. Then, in the L^-norm, S'^ 
produced by the wave-front tracking algorithm is a Cauchy sequence. Denote this unique limit by S 
such that, for any x > 0, §x{U) = lim^_^o Then the map S : [0, oo) xT> ^T> is a uniformly 

Lipschtiz semigroup in L^. In particular, the entropy solution to the initial-boundary problem (jl.ip 
and (|1.7p - ()1.8p constructed by the wave-front tracking algorithm is unique and stable. 

Based on the essential estimates in Sections 3-5, the proof of Theorem 6.1 can be shown in the 
same way as the argument given in [6]. Also see Chen-Li [9]. 

7. Uniqueness of entropy solutions in the viscosity class 

In this section, as an immediate consequence of the estimates obtained in Sections 4-6, we find 
that the semigroup S produced by the wave-front tracking method is the only standard Riemann 
semigroup (SRS) in the sense of Definition 7.1 given below. In other words, the semigroup defined 
by the wave-front tracking method is the canonical trajectory of the standard Riemann semigroup 
(SRS). This yields the uniqueness of entropy solutions in a broader class of viscosity solutions as 
introduced by Bressan in [1]. Furthermore, it coincides with the semigroup trajectory generated 
by the wave-front tracking method. 
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Definition 7.1. Problem (jl.ip and (jl.7p - (|1.8p is said to have a standard Riemann semigroup 
(SRS) if, for some small uq, there exist a continuous mapping 31 : [0,oo) x I? i— )• P and a constant 
L satisfying the following properties: 

(i) (Semigroup property): JioU =U, Jl^j^Jlx^U = 3lr,^+^^U; 

(ii) (Lipschitz continuity): II^^^^C/ — IR^.y||^i < L||C/ — l^H^^i; 

(iii) (Consistency with the Riemann solver): Given piecewise constant initial data U gV, then, 
for all X E [0, vq], the function U(x, ■) = §xU coincides with the solution of p.ip and (jl.7p - 
(|1.8p obtained by piecing together the standard Riemann solutions and the lateral Riemann 
solutions. 

Following the argument in [3], we employ the estimates obtained in Sections 4-6 to conclude 

Theorem 7.1. Suppose that problem (jl.ip and ()1.7p - (|1.8p has a standard Riemann semigroup 
: [0, oo) X D D. Consider the semigroup S produced by the wave-front tracking method, that 
is, §x{U) = lim^^Q §^(U) . Assume C/ G P. Then, for all x > 0, OlxU = SxU. Furthermore, a 
continuous map U : [0, X] D is a viscosity solution of problem (jl.ip and (|1.7p - (|1.8p defined in 
[1] if and only if 

U(x, ■) = JixU for any x G [0,r]. (7.1) 
In particular, a continuous map U : [0, X] D is a viscosity solution if and only if 

U(x,-)=§xU for any X £ [0,T]. (7.2) 

The proof here follows a similar argument to the one presented in [3]. The only difference is 
that there is a strong vortex sheets /entropy waves in our problem. Nonetheless, one can proceed 
with the proof by considering the convergence of the wave-front tracking method which is shown 
in Section 3. 

Remark 7.1. In the simpler cases of the isentropic or isothermal Euler flow (jl.Sp . as well as the 
potential flow, as far as the L^-stability problem is of concern, we realize the same results as for 
the full Euler system (jl.ip . 
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